Gravitomagnetic Accelerators 
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We study a simple class of time-dependent rotating Ricci-flat cylindrically symmetric spacetime 
manifolds whose geodesies admit gravitomagnetic jets. The helical paths of free test particles in 
these jets up and down parallel to the rotation axis are analogous to those of charged particles in a 
magnetic field. The jets are attractors. The jet speed asymptotically approaches the speed of light. 
In effect, such source-free spacetime regions act as "gravitomagnetic accelerators." 

PACS numbers: 04.20.Cv 
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The purpose of this Letter is to point out that general 
relativity in principle permits the existence of dynamic 
source-free spacetime regions in which the speed of free 
test particles can rapidly increase so as to approach the 
speed of light. 

Mechanisms for the acceleration of particles to ultra- 
relativistic speeds are important for the explanation of 
high-energy astrophysical phenomena [TJ [2]. Most of 
the known acceleration processes involve charged parti- 
cles and are electromagnetic in origin [3 . Gravitational 
mechanisms involving arbitrary test particles have been 
the subject of recent investigations [4--8J. We adopt a 
different approach here and present a simple explicit so- 
lution of Einstein's vacuum field equations for which the 
geodesic equation has exact analytic solutions that repre- 
sent gravitomagnetic jets. The jets are attractors and the 
speed of a free test particle in such a jet asymptotically 
approaches the speed of light. 

Consider a warped product spacetime with a metric of 
the form 



ds 2 = e t dZ 2 + e~ t dz 2 , 
where dY? is a 3D stationary metric given by 



dS 2 = ^(_ X 2 d< 2 



1 dr 2 ) + -(-Xdt + d<p) 2 
r 



,-.?, 



(1) 



(2) 



X r- 

and we use units such that c = 1 throughout. Here X r — 
dX/dr and X is the solution of the differential equation 

r 2 X 2 X rr + X r = 0. (3) 

The spacetime metric (fTl) has signature +2 and is Ricci- 
flat, namely, R^ v — [9]. The cylindrical coordinates 
x^ = (t, r, cf>, z) are dimensionless. To transform them 
to physical coordinates x^ — (X't, A _1 r, cf>, Xz), we need 
arbitrary lengthscales A' and A; then, X — X X I 2 X and 
s = As. The conformal factor exp(t), the warping func- 
tion exp(— t) and the gravitomagnetic potential X lead 
to interesting aspects of the gravitational field under con- 
sideration here. 

The cylindrical coordinates in the spacetime domain 
of interest must be admissible. It follows from the Lich- 
nerowicz conditions 110) that the principal minors of the 



metric tensor and its inverse must be negative. Hence, 
XX r > and XQ > 0, where Q = rX r - X. The 
first of these conditions is satisfied if X 2 is a monotoni- 
cally increasing function of r. It turns out that Eq. ([3]) 
is a special case of the generalized Emden-Fowler equa- 
tion. Of the two known explicit solutions of this equa- 
tion, X — constant is not acceptable as metric (fTl) would 
then degenerate into a 3D metric and X = ±(3r/2) _1 / 2 
is such that X 2 is monotonically decreasing with r. In 
fact, with r > taken as radial coordinate, the solutions 
of Eq. pi other than X = constant fall into two classes: 
either X 2 is monotonically increasing or decreasing. In 
this work we focus on the former class, since the latter 
class involves rotating gravitational waves and was inves- 
tigated in detail in Ref. [5]. 

The solutions of Eq. (|3| with X 2 monotonically in- 
creasing all have the general form depicted in Figure [TJ 
If X is a solution of Eq. ^, then so is —X\ henceforth, 
we consider only the X > branch. It is possible to 
show that there is an open set of initial conditions cor- 
responding to admissible solutions X; each exists on a 
radial interval (r^oo), where r^ > with limiting val- 
ues X(rb) = and X r (rb) — oo such that XX r = r^ 2 
at rb- Furthermore, X 2 increases, X 2 decreases and X 2 
approaches infinity as r — > oo. Indeed, r = oo at the 
axis of cylindrical symmetry and r = r^ at the exterior 
boundary cylinder. The physical region of interest S is 
the open hollow expanding cylindrical domain with inner 
boundary around the axis r — oo and outer boundary at 
the null hypersurface r = rv The cylindrical coordinates 
are not admissible at these boundaries, since g tt vanishes 
at r = oo and (— g)~ 1 ' 2 vanishes at r = Tb- 

With admissible cylindrical coordinates in S, we have 
an algebraically general Petrov type I solution of Ein- 
stein's source-free equations with two commuting space- 
like Killing vector fields d^, and d z associated with cylin- 
drical symmetry. The two-parameter isometry group is 
not orthogonally transitive. While d z is hypersurface or- 
thogonal, dej, is not; moreover, the other two coordinates 
(t, r) can be invariantly defined via the magnitudes of 
these Killing vector fields. The t = constant hypersur- 
faces are spacelike, while r = constant hypersurfaces are 



timelike in the physical domain (r(,,oo). Near the axis 
(r —> oo), 



X{r) ~ ar — a 



6ar 2 



6a 2 r 






■oik 



(4) 



where a > and a > are constants. Near the boundary 
(r-»rjj-), 



X(r) 



2 l/2 
^6 



(r-r b ) 1 / 2 +^(r-r b ) + 0((r-r fc ) 3 / 2 ), (5) 



where „4 is a constant. We note that solutions of Ein- 
stein's gravitational field equations with cylindrical sym- 
metry have been the subject of numerous investigations 
(see [TTHT5] and the references cited therein). 

The proper radial distance from the axis to an event 
with r = ro is given by 



p*/2 



X r \ 1/2 



r 3 X 



dr. 



(6) 



At a finite instant of time t, the proper radial distance 
from the axis to the boundary is finite. It is possible 
to show that the condition of elementary flatness is not 
satisfied near the axis [H]. Within the spacetime region 
of interest S, however, there are no curvature singulari- 
ties. To illustrate this fact, we note that there are four 
algebraically independent scalar polynomial curvature in- 
variants in a Ricci-flat spacetime that can be represented 

as na 



Ji = R f>t w P a - iR r>*p,vpo 



7- rt upaa/3 r> P-v . ■ r> r>paal3 p* M" 

-^-2 — I *'p 1 vpfT r *' Ix aj5 ' iIx pi'p<y JX rt a {3 



These turn out to be real in this case and can be ex- 
pressed as 



-2/ 



1-2 = - 



rX 4 X 2 
3e~ 34 



4rX 5 X? 



{F 2 -rX 2 {F + A)), (9) 

F(F + 2rX 2 ), (10) 



where F — r 2 XX r — 1. These invariants are well behaved 
in the interior of the physical region (rf,,oo) and have 
proper limits at the boundaries. Let us note that Xi and 
I2 both diverge in the infinite past (t — > —00), a situation 
that is consistent with the emergence of the universe from 
a singular state as in the standard cosmological models. 
The spacetime domain S rotates about the axis of 
cylindrical symmetry. To elucidate the gravitomagnetic 
aspects of S, imagine the class of fundamental observers 
in S. They are spatially at rest by definition, with spatial 
frames that we can choose to be along the natural direc- 
tions of the cylindrical coordinate system. A unit gyro 
carried by these observers precesses about the z axis with 
frequency (2XQ)~ l , indicating the presence of a gravit- 
omagnetic field parallel to the z axis [T7]. There is ex- 
perimental evidence for gravitomagnetism; indeed, GP-B 




FIG. 1. Left panel: Plot of the solution doX versus r for 
the differential Eq. Q3p. The initial data are X(l) — 1/2 and 
-Xr(l) = 2. The left endpoint occurs at rb ~ 0.88. This is 
an admissible solution of Eq. pj that allows a jet at rj — 1. 
Right panel: Starting from rest at ro = 0.94, halfway between 
ri, and rj, the geodesic equation is integrated for a free test 
particle as in Fig. [2] For < C<f, < 1.65, the particle is 
attracted to the jet at rj (right arrow), while for Cf, > 1.7 
the particle moves to r;, and exits S (left arrow). Similarly, 
for —5 < Cf, < the particle is attracted to the jet, while for 
C<*> < —6 the particle exits via tv 



has recently measured the exterior gravitomagnetic field 
of the Earth [15] . 

We now turn to the motion of free test particles in this 
gravitational field. Let u M = dx^ jdr be the four-velocity 
vector of the test particle, where r is the proper time 
along its world line. The components of u^ along the 
Killing vector fields will be constants of geodesic motion. 
Thus C z = g za u a is the constant linear momentum (per 
unit mass) of the test particle parallel to the z axis, while 
C<f, — g ( j >a u a is the constant angular momentum (per unit 
mass) of the test particle about the z axis. That is, 



(7) 
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= -1 that dt/dr = V/X, 
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:U l 
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(12) 



Here U = dr/dr and we have assumed that t monotoni- 
cally increases with r along the geodesic world line. Thus 
the geodesic equation reduces in this way to the radial 
equation given by 



dU 
~dr~ 



= -i> 



(13) 



In practice, Eq. ^ must also be solved simultaneously 
with respect to the particle's proper time. 

In analogy with the motion of a charged particle in a 
magnetic field, we look for geodesies that are confined to 
a cylinder of fixed radius r 6 (r&, 00). Thus we set U = 
and dU/dr = 0. Equations (|Ll"|) and (12) yield 



dt _ V 

# ~ X(V + C (t> re~ t )' 



while it follows from Eq. ([13]) that 

dt 1 

d4 



Q±{x r /xy/ 2 ' 



(14) 



(15) 




orthonormal tetrad A M /n such that in (t, r, 0, z) coordi- 
nates 



FIG. 2. The result of integration of the geodesic equation for 
tinrelike geodesies attracted to the jet associated with X given 
in Fig. [T] The initial data at tq — are to = 0, ro = 0.94, 
U = 0, 4>o = 0, z = 0, X(r ) w 0.36 and X r (r ) « 2.84. 
In the left-hand plot the parameters are C^ = — 5 and C z = 
1. The right-hand plot is for Cj, = 1.65 and C z = 1. The 
coordinates are (x,y,z), where x = pcos(p, y = psin<ji and 
p = r -1 ' 2 ; moreover, £ = ln| lnz|. We use z instead of z for 
the sake of clarity. 



These are compatible, the latter with the lower sign only, 
if and only if C^ = and there is a radial coordinate 
r J £ (r&,oo) such that r 2 XX r — 1 for r — rj. In this 
case, the geodesic equation can be solved explicitly for 
a special three-parameter set of geodesies parameterized 
by the initial values of (£, <f>, z) such that r remains con- 
stant at rj . The geodesies in this class exhibit helical 
motions up (C z > 0) and down (C z < 0) on the cylinder 
r = r j parallel to the z axis except for the measure zero 
subset (C z — 0) whose motions remain bounded on cir- 
cular orbits about the z axis. It can be shown that for a 
given X, either there is a unique rj, F(rj) = 0, or there 
is none; moreover, there is a non-empty open subset of 
admissible solutions of Eq. p| , corresponding to A < 
in Eq. (pi), that allow these special geodesies. For a so- 
lution of Eq. pi in this class, we call the set of special 
geodesies a gravitomagnetic jet. 

Numerical experinrents reveal that gravitomagnetic 
jets are attractors, see Figs, [l] and [2j That is, the union 
of special geodesies with a given X and rj is a non- 
compact connected invariant manifold that attracts all 
nearby geodesies. Figure [2] highlights the helical motion 
of a gravitomagnetic jet. We note that these helical mo- 
tions up and down within a double-jet configuration have 
the same orientation; the helical sense is positive in our 
case due to our choice of the X > branch. For recent 
studies of helical motions in astrophysical jets, see [!§] . 

Null geodesies can be treated similarly in S and it turns 
out that they have special helical solutions confined to 
r = r j just as in the case of timelike geodesies. Moreover, 
for i -> oo, special timelike geodesies go over to special 
null geodesies as a simple consequence of the dynamical 
equations of motion. 

To explore this important aspect of gravitomagnetic 
jets further, we consider the speed of jets with respect to 
the fundamental observers. These are endowed with an 



A" 



(t) 



((-gtt)- 1/2 , 0,0,0) 



(16) 



and we can choose the spatial frame A'V-, to be along 
the standard cylindrical coordinate axes. Then u^ L — 
u ^(a)' wnere u ^ = 7(h v ) i s the jet four-velocity 
measured by the fundamental observers. Hence, 7 = 
— MpA^cjs and a straightforward calculation reveals that 



7 = 7o (l + C7 z V) 1/2 . 



(17) 



Here 70 is the Lorentz factor corresponding to circular 
motion with C z = 0; that is, 70 = (1 — /Sq)" 1 / 2 , where 
/3q = r,jX 2 (rj). For C z 7^ and t — > 00, 7 diverges 
exponentially with time; that is, the jet speed rapidly 
approaches the speed of light. The gravitational influ- 
ence of the test particle on the spacetime geometry has 
been neglected in our work; clearly, this approximation 
eventually breaks down for the gravitomagnetic jets. 

Equation pi) is invariant under the scale transforma- 
tion (r,X) n- (f,X), where f = or and X = o~ 1 / 2 X 
for a G (0,oo). This scale invariance can be used to 
reduce Eq. pi) to a first order system of autonomous 
equations. For instance, let 6 = — hxX and consider 
scale-invariant variables p and q, where p(9) = X/(rX r ) 
and q(9) = (F + l)^ 1 ; then, Eq. ^ reduces to the Lotka- 
Volterra system 

!=K,-*-D, § = «(2p-« + l). (18) 

A solution X is admissible when < p < 1 for all 9; 
similarly, an admissible solution allows jets when q(9) = 1 
for some 9 > 0, where 9 = —00 at the symmetry axis and 
9 = 00 at the outer boundary. A detailed investigation 
of system ( 18 ) — for admissible solutions that allow jets — 
reveals that p/q = rX 2 for q = 1 lies in the interval 
(0,j), where j ps 0.4. This means that ^0, the speed of 
the free test particle on a circle of radius rj, is such that 
A) G (0, j 1/2 ), where j 1 / 2 w 0.63. 

It follows from these results that free test particles in a 
gravitomagnetic jet can in principle start out with speeds 
near zero, but they then inevitably undergo rapid "accel- 
eration" to almost the speed of light. Similarly, one can 
study the motion of other geodesies in S with respect to 
the fundamental observers. In the simple numerical ex- 
periments whose results are presented in Figs.[l]and[2] we 
followed the variation of the Lorentz factor from the ini- 
tial point to late times. There is indeed a vast difference 
between the geodesies that leave S via r& and those that 
are attracted to the jet — see the right panel of Fig. [I] In 
the former case, the Lorentz factor initially decreases but 
then increases as the geodesic exists 5, remaining within 
about an order of magnitude of unity. In the latter case, 
the Lorentz factor quickly diverges to infinity — starting 
from sa 3.5 (for the left-hand plot) and w 2.9 (for the 
right-hand plot) in Fig. [2] 



The existence of a jet is a scale-invariant property; that 
is, for every admissible solution with a jet, scaling leads 
to a one-parameter family of solutions of the same kind, 
where the jet now occurs at orj. One can use this prop- 
erty to set rj = 1 for every jet solution. Thus Eq. (J3J) 
can be integrated with initial conditions that at Tj = 1, 
X(l) = i?" 1 and X r (l) = tf; in this case, r b « 1 - t?~ 2 /2 
for i? ^ 1. The solution is admissible once /3 2 = $ -2 < j; 
that is, $ > j -1 / 2 ~ 1.6. For example, in Fig. [l] we have 
■0 = 2 and hence /3o = 0.5 for the corresponding jet de- 
picted in Fig. [2] A detailed treatment of gravitomagnetic 
jets is contained in Ref. [2D] . 

The empty annular region of physical interest that al- 
lows gravitomagnetic jets could in principle be joined to 



two concentric circular cylindrical regions with matter 
such that the inner one surrounds the axis (r = oo) and 
the outer one has its inner boundary at r = rv Our 
treatment has been greatly simplified by the assumption 
of cylindrical symmetry; nevertheless, we hope that in 
the context of general relativity similar gravitomagnetic 
accelerators may emerge under physically more realistic 
circumstances. In any case, this appears to be a promis- 
ing approach in the search for the origin of high-energy 
astrophysical jets. Further recent results in this direction 
are contained in Ref. [2"T] . 
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